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Abstract. We investigate cosmological predictions on the early uni- 
verse based on the noncommutative geometry models of gravity coupled 
to matter. Using the renormalization group analysis for the Standard 
Model with right handed neutrinos and Majorana mass terms, which 
is the particle physics content of the most recent noncommutative ge- 
ometry models, we analyze the behavior of the coefficients of the grav- 
itational and cosmological terms in the Lagrangian derived from the 
asymptotic expansion of the spectral action functional of noncommuta- 
tive geometry. We find emergent Hoyle-Narlikar and conformal grav- 
ity at the see-saw scales and a running effective gravitational constant, 
which affects the propagation of gravitational waves and the evapora- 
tion law of primordial black holes and provides Linde models of negative 
gravity in the early universe. The same renormalization group analysis 
also governs the running of the effective cosmological constant of the 
model. The model also provides a Higgs based slow-roll inflationary 
mechanism, for which one can explicitly compute the slow-roll parame- 
ters. The particle physics content allows for dark matter models based 
on sterile neutrinos with Majorana mass terms. 
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1. Introduction 



The idea of using noncommutative geometry to give a conceptual math- 
ematical formulation of the standard model of elementary particle physics 
dates back to the work of Connes |l9j . It was shown more recently in [T7] 
(see also [20], [JJ, and Chapter 1 of [22]) that the noncommutative geometry 
model of particle physics can be made compatible with right handed neutri- 
nos and neutrino masses and that the full Lagrangian of the standard model 
with Majorana mass terms for the right handed neutrinos can be derived by 
a computation from a very simple input of an almost commutative space, 
the product of an ordinary spacetime manifold and a finite noncommutative 
space. The resulting physical Lagrangian in this noncommutative geometry 
model is obtained from the asymptotic expansion in the energy scale A of 
a natural action functionals defined on noncommutative spaces, the spectral 
action, [2] . Among the most interesting features of these models of particle 
physics based on noncommutative geometry is the fact that the physical 
Lagrangian of the model is completely computed from a simple geometric 
input (the choice of a finite dimensional algebra) , so that the physics is very 
tightly constrained by the underlying geometry. For reasons of space, we 
cannot include here any introductory material about Noncommutative Ge- 
ometry, but we suggest the interested readers to look at the survey paper 
|23] for a user-friendly introduction based on examples, as well as the books 
|21] and [22] for a more complete treatment. 

We focus here on the noncommutative geometry model obtained in |17j . 
The corresponding physical Lagrangian computed from the asymptotic ex- 
pansion of the spectral action contains the full Standard Model Lagrangian 
with additional Majorana mass terms for right handed neutrinos, as well as 
gravitational and cosmological terms coupled to matter. The presence of 
these terms and their relation to the particle physics content of the model 
make this approach of interest to theoretical cosmology. The gravitational 
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terms include the Einstein-Hilbert action with a cosmological term, a topo- 
logical term related to the Euler characteristic of the spacetime manifold, 
and, additionally, a conformal gravity term with the Weyl curvature tensor 
and a conformal coupling of the Higgs field to gravity. Without these last 
two contributions essentially one would be dealing with the usual minimal 
coupling of the standard model to gravity, but the presence of an additional 
non-minimal conformal coupling has relevance to various cosmological mod- 
els that have been the object of recent investigation, especially in the context 
of modified theories of gravity. Another important way in which this model 
differs from the minimal coupling of gravity and matter is the dependence 
of the coefficients of the gravitational terms upon the Yukawa parameters of 
the particle physics content of the model. This feature, which is our main fo- 
cus of investigation in the present paper, is unique to these noncommutative 
geometry models and does not have an analog in other particle physics and 
cosmology models obtained from geometric settings such as string theories, 
extra dimensions, or brane worlds. Some conceptual similarities with these 
other approaches exist though, in the sense that in noncommutative geom- 
etry models one typically modifies ordinary spacetime by taking a product 
with a noncommutative space and this may be thought of as another possible 
way to enrich it with extra dimensions. 

The fact that, as shown in |17j . the model lives naturally at unification 
scale, means that in cosmological terms it provides us with early universe 
models, hence it is interesting in terms of possible inflationary mechanisms. 
Extrapolations to lower energies are possible using renormalization group 
analysis, though extensions to cosmological models of the more recent uni- 
verse only become possible when non-perturbative effects in the spectral 
action are also taken into account. The main motivation for considering 
these noncommutative geometry models in a cosmological context is that 
the nontrivial dependence of the cosmological and gravitational parameters 
on the particle physics content is, as we mentioned above, significantly differ- 
ent from other physical models, hence likely to provide inflationary scenarios 
in the early universe that differ significantly from other models. 

This paper is the first of a planned series dedicated to an investigation 
of the cosmological implication of the noncommutative geometry models in 
particle physics. In the present paper, we concentrate on a renormalization 
group analysis of the coefficients of the gravitational terms in the action. In 
fact, the asymptotic formula for the spectral action used in [T7] shows that 
these coefficients are functions of certain parameters, which in turn depend 
on the data of the Yukawa parameters of the standard model with Majorana 
mass terms for right handed neutrinos. They also depend on three additional 
parameters of the model, one of which is fixed by a grand unification type 
condition on the coupling constants. 

The Yukawa parameters run with the renormalization group equations 
of the particle physics model. In particular, since the noncommutative ge- 
ometry model lives naturally at unification scale, one can input boundary 
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conditions at that energy scale and follow the renormalization group equa- 
tions towards lower energies and investigate the effect of this running on the 
gravitational part of the model. One expects that, when running towards 
lower energies, nonperturbative effects in the spectral action will progres- 
sively become non-negligible. This can limit the range of validity of this type 
of argument based on the asymptotic expansion alone, and on renormaliza- 
tion group analysis to cosmological models for the very early universe, that 
is, for sufficiently high energies where the asymptotic expansion holds. Any 
extrapolation to the modern universe would then have to take into account 
the full spectral action and not only its asymptotic form. 

In the present paper we focus on early universe models and on the as- 
ymptotic form of the spectral action. For the renormalization analysis, we 
rely on a detailed study of the renormalization group equations (RGE) for 
the extension of the standard model with right handed neutrinos and Ma- 
jorana mass terms carried out in [I], even though their choice of boundary 
conditions at unification is different from some of the boundary conditions 
assumed in [T7]. The boundary conditions proposed in [1] are dictated by 
particle physics considerations, while some of the constraints considered in 
|17] came from analyzing particular geometries, such as the flat space case. 
For example, for simplicity the Majorana masses were assumed in [17J to be 
degenerate, all of them close to unification scale, while here we are mostly 
interested in the non-degenerate case, with three different see-saw scales 
between the electroweak and unification scales, which leads to a more inter- 
esting behavior of the gravitational terms in the model. We plan to return 
to a more general analysis of the RGE flow of [1] with a wider range of 
possible boundary conditions in followup work. 

The RGE flow of [1] runs between a unification energy, taken there to be 
of the order of 2 x 10 16 GeV, down to the electroweak scale of 100 GeV. In 
terms of cosmological timeline, we are looking at the behavior of the model 
between the unification and the electroweak era. This means that, in terms 
of matter content, only the Higgs field and its coupling to gravity is relevant, 
so we mostly concentrate on the part of the Lagrangian of [T7] that consists 
of these terms. Since this era of the early universe is believed to include 
the inflationary epoch, we look especially at different possible inflationary 
scenarios provided by this noncommutative geometry model. 

Our main results in this first paper are to show that, using the information 
on the particle physics content, it is possible to obtain cosmological models 
of the early universe with variable gravitational and cosmological constant, 
hence providing a range of different mechanisms for inflation, realized by 
the running of the effective gravitational constant and by its coupling to 
the Higgs field, or by the running of the effective cosmological constant of 
the model, or by a combination of these. We also show phenomena where, 
near particular energy scales and for special geometries, the usual Einstein- 
Hilbert action ceases to be the dominant contribution and the model comes 
to be dominated, at certain scales, by conformal gravity and an emergent 
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Hoyle-Narlikar cosmology. We discuss how the running of the gravitational 
parameters of the model influences the behavior of the evaporation of pri- 
mordial black holes by Hawking radiation. While the type of effects that we 
see in this model, which depend on the presence of variable effective grav- 
itational and cosmological constants, are qualitatively similar to scenarios 
of negative gravity in the early universe previously analyzed in theoretical 
cosmology (0, @, p], [Mj, EE [26], |35], [36], [37], [39], 00], gl], 03], 
|48] ) . the mechanism that produces these effects in the noncommutative ge- 
ometry model is substantially different from those described in these earlier 
references, which makes the quantitative behavior also different and distin- 
guishable from other models. In fact, most of the effects we investigate in 
this paper depend directly on the expression of the coefficients of the gravita- 
tional and bosonic terms in the asymptotic expansion of the spectral action 
in terms of the Yukawa parameters of the underlying particle physics model. 
This is a purely geometric property of this model and it comes directly from 
the presence of the "small extra dimensions" in the form of the zero dimen- 
sional (but K-theoretically six dimensional) finite noncommutative space in 
addition to the extended spacetime dimensions. 

While the energy range where the renormalization group analysis applies 
limits the results based only on the perturbative expansion of the spectral 
action to early universe models, if some of the results obtained in this paper 
persist when non-perturbative effects in the spectral action become signif- 
icant, they may provide possible dark energy and dark matter predictions. 
For instance, the behavior of the variable effective cosmological constant 
may lead to dark energy scenarios in the more recent universe. Moreover, 
we show that the particle physics content of the model is consistent with 
dark matter models based on right handed neutrinos with Majorana mass 
terms in |33j, [S], [H]. In fact, the particle content is the same as in the 
z^MSM model with three active and three sterile neutrinos. What is needed 
in order to relate the model of [17^ to these dark matter models is a choice of 
boundary conditions that makes it possible for at least one of the Majorana 
masses to descent to somewhere near the electroweak scale, hence providing 
sterile neutrinos with the characteristics required to give acceptable dark 
matter candidates (see the analysis in [33 ) . We will return to a closer anal- 
ysis of dark energy and dark matter implications of the noncommutative 
geometry models in a planned continuation of this work. 

In a related but different direction, recent work on some cosmological 
aspects of the model of [T7] was also done in [38J . 
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2. The asymptotic formula for the spectral action and the 
gravitational parameters 

We recall here briefly the main aspects of the noncommutative geometry 
model of particle physics derived in [T7] that we need to use in the rest of 
the paper. We refer the reader to [17 and to Chapter 1 of [22 for a detailed 
treatment. The reader who wishes to skip this preliminary part can start 
directly with the asymptotic expansion of the spectral action recalled in 
§2.5[ which is what we concentrate on in the rest of the paper, but we prefer 
to add a few words on the derivation of the model via noncommutative 
geometry for the sake of completeness. 

2.1. Spectral triples and the spectral action functional. The particle 
physics models based on noncommutative geometry, both the original one 
of [19] and the new one of [T7] that incorporates right handed neutrinos and 
neutrino mixing with Majorana mass terms, are based on the formalism of 
spectral triples. These were introduced by Connes [18] as an extension of 
the notion of Riemannian manifold to noncommutative geometry. The data 
(A,TL,D) defining a (real) spectral triple are summarized as follows. 

• A is an involutive algebra with unit. Requiring the algebra to be 
unital corresponds to working with compact manifolds. (Extensions 
of the notion of spectral triple to non-unital cases have also been 
developed.) 

• TL is a separable Hilbert space endowed with a representation tt : 
A — ► £>(H) of the algebra A by bounded linear operators. 

• D = D' is a linear self-adjoint operator acting on TL. Except for 
finite dimensional cases, D is in general not a bounded operator, 
hence it is only defined on a dense domain. 

• D has compact resolvent: (1 + D 2 )~ l l 2 is a compact operator. 

• The commutators [Tr(a),D] are bounded operators for all a £ A. 

• The spectral triple is even if there is on If a Z/2-grading 7 satisfying 
[7, n(a)] = and Dj = — 7.D. 

• The spectral triple has a real structure if there is an antilinear iso- 
morphism J : TL — > TL with J 2 = e, JD = e'DJ, and J7 = e"^/J, 
where the signs e, e', and e" determine the KO-dimension modulo 8 
of the spectral triple, according to by the table 



n 





1 


2 


3 


4 


5 


6 


7 


e 


1 


1 


-1 


-1 


-1 


-1 


1 


1 


e' 


1 


-1 


1 


1 


1 


-1 


1 


1 


e" 


1 




-1 




1 




-1 





• The Hilbert space TL has an ,4-bimodule structure with respect to the 
action of A defined by b° = Jb* J -1 and satisfying the commutation 
condition [a, b°] = for all a and b in A. 

• The operator D satisfies the order one condition [[D,a],6°] = 0, for 
all a, b 6 A- 
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Commutative geometries, which in this context means ordinary Riemann- 
ian manifolds, can be described as spectral triples: for a compact spin Rie- 
mannian manifold X the associated spectral triple (C°°(X), L 2 (X, S), Dx) 
is given by the algebra of smooth functions, the Hilbert space of square inte- 
grable spinors, and the Dirac operator. The metric tensor can be recovered 
from these data. For an even dimensional manifold "fx = 75 is the grading 
H = 7i + © 7i~ on the spinor bundle given by the usual chirality operator 75 
and the real structure Jx is the charge conjugation operator. Examples of 
spectral triples associated to objects that are not manifolds include a wide 
range of geometries such as quantum groups, fractals, or noncommutative 



tori. As we recall in {2.2 below, the spectral triples involved in the parti- 
cle physics models are of a very special form which is almost commutative, 
namely a product of an ordinary manifold with a small noncommutative 
space. 

It was shown by Chamseddine and Connes |14| that there is a natural 
action functional on a spectral triple. This spectral action functional is 
defined as Tr(/(Z)/A)), where / > is a cutoff function and A is the energy 
scale. There is an asymptotic formula for the spectral action, for large energy 
A, of the form 

(2.1) Tr(/(L>/A))~ Yl fkA k f\D\- k + f(0)CD(0) + o(l), 

fcSDimSp 

where = J °° f{v)v k ~ l dv are the momenta of the function / and the 
noncommutative integration is defined in terms of residues of zeta functions 

(2.2) £ a>D (s) = ^(a\D\- s ). 



The sum in (2.1) is over points in the dimension spectrum of the spectral 



triple, which is a refined notion of dimension for noncommutative spaces, 



consisting of the set of poles of the zeta functions (2.2 ). 



2.2. The noncommutative space of the model. The main result of 
|17] is a complete derivation of the full standard model Lagrangian with 
additional right handed neutrino, lepton mixing matrix and Majorana mass 
terms, by a calculation starting from a very simple geometric input. The 
initial ansatz used in [T7] is the choice of a finite dimensional algebra, the 
left-right symmetric algebra 

(2.3) A LR = C(BU L (BM R (BM 3 {C), 

where Hi and Mr are two copies of the real algebra of quaternions. 

The representation is then naturally determined by taking the sum Ai of 
all the inequivalent irreducible odd spin representations of Alr, so that the 
only further input that one needs to specify is the number N of generations. 
The (finite dimensional) Hilbert space is then given by N copies of A4 , 

H F = (B N M. 
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The Hilbert space with the Alr action splits as a sum Tip = Hf © Hf of 
matter and antimatter sectors, and an orthogonal basis olTif gives all the 
fermions of the particle physics model 



(2.4) 



V L = 




<g>i° 


vr = 




®l° 


ZL = 




8)1° 


ZR = 


li>fl 




UL = 


IT)l 


®3° 


UR = 


1 T>« 


® 3° 


d L = 


1 ih 


® 3° 


d,R = 


1 Dr 


3° 



respectively giving the neutrinos, the charged leptons, the u/c/t quarks, and 
the d/s/b quarks in terms of the representation of Alr- Here | |) and | J.) 
are the basis of the 2 representation of H where the action of A G C C H 
is, respectively, by A or A, and 1° and 3° are the actions of C and M$(C), 
respectively, through the representation a = Ja* J~ l . 

The Z/2-grading jp exchanges the left and right chirality of fermions 
and the real structure operator Jp exchanges the matter and antimatter 
sectors and performs a complex conjugation. These properties of jf and Jf 
suffice to determine the KO-dimension modulo 8 of the resulting spectral 
triple and an interesting aspect is that, unlike in the earlier particle physics 
models based on noncommutative geometry, in this case the KO-dimension 
is 6 modulo 8, although the metric dimension is zero. 

The order one condition on the Dirac operator is seen in |17j as a coupled 
equation for a subalgebra _4p C Alr and a Dirac operator and it is shown 
that there is a unique subalgebra of maximal dimension that allows for the 
order one condition to be satisfied. The algebra Af is of the form 

(2.5) A F = C©M©M 3 (C), 

where the first summand embeds diagonally into C©H in Alr, thus breaking 
the left-right symmetry. It is expected, though presently not known, that 
this symmetry breaking should be dynamical. This geometric argument 
identifying the maximal algebra on which the order one condition can be 
satisfied was later extended in [15] to more general ansatz algebras than 
Alr, but with the same resulting Af- 

2.3. Dirac operators: Yukawa parameters and Majorana masses. 

The selection of the subalgebra .4.^ for the order one condition for the Dirac 
operator is what produces geometrically in this model the Majorana mass 
terms for right handed neutrinos. In fact, one has in [T7] a complete classifi- 
cation of the possible Dirac operators on the noncommutative space (Af, H) 
compatible with jf and Jf (see also |10] for a more general discussion of 
moduli spaces of Dirac operators for finite spectral triples). These are all of 
the form 
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with S = Si © (S3 © I3) and T = Yr : \ur) — > Jp \vr), and with Si and S3 
respectively of the form 



( 





r (Ti) 





\ 











Y ] 




y (Ti) 













V 


%1) 








/ 


( ° 





r (T3) 





\ 











Y ] 




y (T3) 













V 










/ 



Here the N x iV-matriees involved in the expression of Si and S3 are the 
Yukawa matrices that give Dirac masses and mixing angles. These are ma- 
trices in GL3(C) in the case oi N = 2> generations: Y e = Yn\\ is the Yukawa 
matrix for the charged leptons, Y v = Vm) for the neutrinos, Y& = Yng\ 
for the d/s/b quarks, and Y u = Yn$\ for the u/c/t quarks. Moreover, the 
remaining term M = Yj£, with T denoting transposition, gives the matrix T 
in D(Y) and is the symmetric matrix of the Majorana mass terms for right 
handed neutrinos. 

Thus, the model of [T7j has three active and three sterile neutrinos as in 
the vMSM model, see [33], @3], @5], though in |T7], unlike in the vMSM 
model, it is assumed that the three sterile neutrinos all have masses well 
above the electroweak scale. The see-saw relation Yj M~ 1 Y u for neutrino 
masses is obtained in [T7] geometrically from the fact that the restriction of 
the Dirac operator D(Y) to the subspace of Jip spanned by ur, i>l, z/r, vl 
is of the form 

/ Ml M r \ 

(2.6) M » ° ° ° f , 

y ' M R Ml 

V My / 

where M u is the neutrino mass matrix, see Lemma 1.225 of [22] . We return 
to discuss the relation of the model of [T7] to the z/MSM model of |H], [15] 
and to other sterile neutrinos scenarios of [33 in the context of dark matter 
models in cosmology, see §[5] below. 

The spectral triple that determines the physical Lagrangian of the model 
through the asymptotic expansion of the spectral action is then the product 
geometry X x F, of a 4-dimensional spacetime X, identified with the spec- 
tral triple (C°°(X), L 2 (X, S), Dx), and the finite noncommutative space F 
defined by the spectral triple (Af,T~If, Dp) with Dp of the form D(Y) as 
above. The product is given by the cup product spectral triple (A,TC,D) 
with sign 7 and real structure J 

• A = C 00 (X)®A F = C 00 (X,A F ) 



10 



MATILDE MARCOLLI AND ELENA PIERPAOLI 



• H = L 2 (X, S)®H F = L 2 (X, S <g> H F ) 

• D = D x (8> 1 + 75 <8> D F 

• J = Jx 8> Jf and 7 = 75 <g> 7^. 

The action functional considered in jTTj to obtain the physical Lagrangian 
has a bosonic and a fermionic part, where the bosonic part is given by the 
spectral action functional with inner fluctuations of the Dirac operator and 
the fermionic part by the pairing of the Dirac operator with fermions, 

(2-7) Tr(f(D A /A)) + ^(JlD A £). 

Here Da = D + A + e' J AJ' 1 is the Dirac operator with inner fluctuations 
given by the gauge potentials of the form A = A* = J2 k ak[D,bk], for ele- 
ments ak,bk G A. The fermionic term (J£, Da£) should be seen as a pairing 
of classical fields £ G 7i + = {£ G Ti \ 7<£ = £}, viewed as Grassman variables. 
This is a common way of treating Majorana spinors via Pfaffians, see §16.2 
of [22]. 



While this fermionic part is very important for the particle physics content 
of the model, as it delivers all the fermionic terms in the Lagrangian of the 
Standard Model, for our purposes related to cosmological models of the early 
universe, it will suffice of concentrate only on the bosonic part of the action, 
given by the spectral action term Tr(/(D^/A)), since during a good part of 
the cosmological period between the unification and the electroweak epoch 
the Higgs field is the matter content that will be mostly of relevance, |30j. 



2.4. Parameters of the model. As we have recalled above, the geomet- 
ric parameters describing the possible choices of Dirac operators on the 
finite noncommutative space F correspond to the Yukawa parameters of the 
particle physics model and the Majorana mass terms for the right handed 
neutrinos. We recall here some expressions of these parameters that appear 
in the asymptotic expansion of the spectral action and that we are going to 
analyze more in detail later in the paper. We define functions a, b, c, d, e of 
the matrices Y u , Y^, Y u , Y e and of the Majorana masses M in the following 
way: 

a = Tt(yJy v + Yhe + KYuYu + Y ] d Y d )) 

b = Tr((yjy,) 2 + (Yhe) 2 + 3(yJy u ) 2 + 3(y]y d ) 2 ) 

(2.8) c= Tr(MMt) 

8= Tr((MMt) 2 ) 

e= Tr(MM^Y v ). 
In addition to these parameters, whose role we describe in §2.5| we see clearly 



from (2.1) that the asymptotic formula for the spectral action depends on 
parameters given by the momenta of the cutoff function / in the spectral 
action. Since the noncommutative space here is of the simple form X x F, 
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the only contributions to the dimension spectrum, hence to the asymptotic 
formula for the spectral action come from three parameters /o, /2, fi, where 
f = /(0) and for k > 

poo 

fk = / f(v)v k ~ l dv. 
Jo 

2.5. The asymptotic expansion of the spectral action. It was proved 
in [T7] that the asymptotic formula (2.1) applied to the action functional 
Tr(f(DA/A.)) of the product geometry X x F gives a Lagrangian of the 
form 

S= ^(48/ 4 A 4 -/ 2 A 2 c + ^9) J J-gd^x 
96/ 2 A 2 -/ c f 4 

+ ^ J R Vad x 

+ / ( ¥ R * R * ~ 3 ^ 

+ ( - L/2A 2 2 + e/0) /M 2 ^ 

(2.9) 



+ / I^V>lW* 



12 7T 2 

2vr 2 



R\<p\ 2 ^d 4 x 

\<p\ 4 Vgd 4 x 



We see from this expansion how the coefficients of all the terms in this 
resulting action functional depend on the Yukawa and Majorana parameters 
through their combinations of the form o, b, c, d, e defined as in (2.8), and 
from the three additional parameters /o, / 2 , fi- 

The term of (2.9) with the Yang-Mills action for the gauge bosons, 

^2 / {9lG^G^+glF« u F^+ 5 -glB, u B^)V9d 4 x, 

contains the coupling constants g\, g 2 , <?3 of the three forces. As shown in 
|17] , the standard normalization of these Yang-Mills terms gives the GUT 
relation between the three coupling constants and fixes the fact that this 
model lives naturally at a preferred energy scale given by the unification 
scale A = A un i t . The normalization of the Yang-Mills terms fixes the value 
of the parameter /o to depend on the common value g at unification of the 
coupling constants: as shown in §4.5 and §5.1 of |17j one obtains 
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One also normalizes the kinetic term for the Higgs as in [T7] by the change 
of variables H = to get \ J \DH\ 2 ^/gd 4 x. 

The normalization of the Yang-Mills terms and of the kinetic term of the 
Higgs then gives, at unification scale, an action functional of the form 



+ a y C^C^y/gd^x + ro J R*R*^d 4 x 

(2.11) + ±f \DH\ 2 ^d 4 x-fi 2 J \H\ 2 ^d A x 

- $oj R\H\ 2 ^d 4 x + X J \H\ 4 ^d 4 x 

+ | J + F% v F>™ + B, u Bn ^ d 4 x, 

where the coefficients are now 

96/ 2 A 2 - f c 



l 

72 



2kI 247T 2 

70= ^(48/ 4 A 4 -M 2 C +^i 

7T Z 4 

3/o 



(2.12) 



,2 



IOtt 2 
ll/o 
60vr 2 

A- ^-l 

Jo a 

12 



6> = 

Ao 



^ 2 b 



2/oa 2 ' 

again as a function of the Yukawa and Majorana parameters through the 
coefficients a, b, c, 0, e of (2.8), and of the two remaining free parameters of 
the model, ji and 74, after the value of fo has been fixed by the unification 
condition. 



3. Renormalization group and running parameters 

All the Yukawa parameters Y u , Y^, Y u , Y e , as well as the Majorana mass 
terms M are subject to running with the renormalization group equations 
(RGE) dictated by the particle physics content of the model, in this case 
the Standard Model with additional right handed neutrinos with Majorana 
mass terms. Consequently, also the parameters a, b, c, t>, and e of (2.8) run 
with the renormalization group flow as functions of A, with assigned initial 
conditions at A = A nn jj, which is the preferential energy scale of the model. 
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Some estimates based on renormalization group analysis were obtained 
already in |17j . for the Higgs and the top quark masses, but those were 
based, in first approximation, on just the renormalization group equations 
at 1-loop for the minimal Standard Model. 

In this section we analyze the running of the parameters of the model 
with the renormalization group flow, using the full RGE of the extension of 
the Standard Model by right handed neutrinos and Majorana masses. There 
is an extensive literature available in particle physics on the relevant RGE 
analysis, see for instance [2], [3], [12]. We use here a more detailed analysis 
of the renormalization group flow, again to one-loop order, for the Standard 
Model with additional Majorana mass terms for right handed neutrinos, 
as given in [I] and implemented by the authors of [I] in the Mathematica 



package http://www.ph.tum.de/~ rge/REAP/. 

The full renormalization group equations for this particle physics model 
have beta functions given by 

(3.1) 16tt 2 P 9i = bi gf with (b su{3) , b su{2) , b u{1) ) = (-7, -y, 

where [I] is using here a different normalization from [T7| and the factor 5/3 
has been now included in g\. Thus, as for the minimal Standard Model, at 
1-loop order the RGE for the coupling constants uncouple from those of the 
other parameters. We then have for the Yukawa matrices 

(3.2) 16vr 2 fa = Y U {\Y^Y U - *Y}Y d + a - ^ - \g\ - 8gl) 

(3.3) 167T 2 (3 Yd = Y^yJy, - \y^Y u + a - \g\ - \g\ - 89%) 

(3.4) 16vr 2 = Y U ( 3 2 Y}Y U - \Y}Y e + a - ^g\ - \&) 

(3.5) 16vr 2 p Ye = Y e (^Y e - 3 2 Y}Y U + a - \g\ - \g\). 

The RGE for the Majorana mass terms has beta function 

(3.6) 16vr 2 M = Y v YlM + M(Y u Yjf 
and the one for the Higgs self coupling A is given by 

(3.7) 167T 2 f3\ = 6A 2 - 3A(3 5 2 + 3 gf) + 3<? 2 4 + ^g\ + 5 2 ) 2 + 4Ao - 8b 

In the treatment of the renormalization group analysis given in the refer- 
ences mentioned above, one assumes that the Majorana mass terms are non- 
degenerate, which means that there are different see-saw scales at decreasing 
energies in between unification and the electroweak scale. In between these 
see-saw scales, one considers different effective field theories, where the heav- 
iest right handed neutrinos are integrated out when one passes below the 
corresponding see-saw scale. In practice the procedure for computing the 
RGE for this type of particle physics models can be summarized as follows. 
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Run the renormalization group flow down from unification energy 
A un i/ to first see-saw scale determined by the largest eigenvalue of 
M, using assigned boundary conditions at unification. 

(3) 

Introduce an effective field theory where Y) is obtained by remov- 
ing the last row of Y v in the basis where M is diagonal and is 
obtained by removing the last row and column. 
Restart the induced renormalization group flow given by the equa- 
tions d3T2b— ([377b with Y u replaced by Y"J 3) and M replaced by 



and with matching boundary conditions at the first see-saw scale. 
Run this renormalization group flow down to second see-saw scale. 

• Introduce a new effective field theory with Y^ and obtained 
by repeating the above procedure starting from Y"J 3 ^ and 

• Run the induced renormalization group flow for these fields with 
matching boundary conditions at the second see-saw scale, down 
until the first and lowest see-saw scale. 

• Introduce again a new effective field theory with YS and and 
matching boundary conditions at the first see saw scale. 

• Run the induced RGE down to the electroweak energy A = A ew . 

The procedure illustrated here assumes that the three see-saw scales are all 
located between the unification and the electroweak scale, that is, that all 
the sterile neutrinos are heavy. 

Notice that there is a difference between the boundary conditions assumed 
in [1] for M at unification energy and the assumption made in [T7] on the 
Majorana mass terms. In fact, in §5.5 of [T7] under the assumptions of 
flat space and with Higgs term \H\ 2 sufficiently small, it is shown that one 
can estimate from the equations of motion of the spectral action that the 
largest Majorana mass can be as high as the unification energy, while in [1] 
the unification scale is taken at around 10 16 GeV but the top see-saw scale 
is around 10 14 GeV. For the purpose of the present paper we only work 
with the boundary conditions of [1 , while a more extensive study of the 
same RGE with a broader range of boundary conditions will be considered 
elsewhere. We report the explicit boundary conditions of [1] at A = A un if 
in the appendix. 

3.1. Running parameters and see-saw scales. In the following, we as- 
sume as in [IT] that the value of the parameter /q is fixed by the relation 



(2.10). In terms of the boundary conditions at unification used in [1], this 



gives as value for either one of the following 

(3.8) f = ^ = 8.52603, f = ^ = 9.46314, f = ^ = 9.36566. 

These three different choices come from the fact that, as it is well known, 
the values for the three coupling constants do not exactly meet in the min- 
imal Standard Model, nor in its variant with right handed neutrinos and 
Majorana mass terms. Notice that here we are already including the factor 
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Figure 2 . The coefficient c as a function of the energy scale 
A near the three see-saw scales. 



of 5/3 in gf, unlike in [T7], where we have g\ = g\ = |<7 2 at unification. 
We shall perform most of our explicit calculations in the following using the 
first value for /o and it is easy to check that replacing it with either one of 
the others will not affect significantly any of the results. The parameters fi 
and remain free parameters in the model and we discuss in ^4] how they 
can be varied so as to obtain different possible cosmological models. 

We concentrate here instead on the coefficients a, b, c, d, and e of (2.8 ) and 
on their dependence on the energy scale A through their dependence on the 
Yukawa parameters and the Majorana mass terms and the renormalization 
group equations (3.2)-(3.7). The renormalization group flow runs between 
the electroweak scale A ew = 10 2 GeV and the unification scale, chosen as in 
[I] at 2 x 10 16 GeV. 

By solving numerically the equations and plotting the running of the 
coefficients (2.8) one finds that the coefficients o and b show clearly the 
effect of the first (highest) see-saw scale, while the effect of the two lower 
see-saw scales is suppressed. 

The running of the coefficients c and Z> exhibits the effect of all three see- 
saw scales, while the running of the remaining coefficient e is the only one 
that exhibits a large jump at the highest see-saw scale. 
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Figure 3 . The coefficient D as a function of the energy scale 
A near the three see-saw scales. 




Figure 4. The coefficient e as a function of the energy scale 
A near the three see-saw scales. 



Notice that the lack of differentiability at all the see-saw scales is in- 
evitable, due to the procedure used in [1] and recalled here above for the 
construction of the effective field theories in between the different see-saw 
scales. 



4. COSMOLOGICAL IMPLICATIONS OF THE MODEL 



We now use the information on the running of the coefficients (2.8) with 
the renormalization group to study the effect on the coefficients of the gravi- 
tational and Higgs terms in the asymptotic expansion of the spectral action. 
We derive some information about cosmological models of the early universe 
that arise naturally in this noncommutative geometry setting. In particular, 
we focus on the following aspects. 



EARLY UNIVERSE MODELS FROM NONCOMMUTATIVE GEOMETRY 



17 



• Spontaneously arising Hoyle-Narlikar cosmologies in Einstein-Hilbert 
backgrounds. 

• Linde's hypothesis of antigravity in the early universe, via running 
of the gravitational constant and conformal coupling to the Higgs. 

• Gravity balls from conformal coupling to the Higgs field. 

• Detectable effects on gravitational waves of the running gravitational 
constant, as in modified gravity theories. 

• Primordial black holes with or without gravitational memory. 

• Higgs based slow-roll inflation. 

• Varying effective cosmological constant and vacuum-decay. 

• Cold dark matter from Majorana masses of right handed neutrinos. 

The main features of the noncommutative geometry model that will be 
discussed in the following and that lead to the effects listed above are sum- 
marized as follows. 

• Variable effective gravitational constant. 

• Variable effective cosmological constant. 

• Conformal gravity. 

• Conformal coupling of the Higgs field to gravity. 

4.1. Einsten gravity and conformal gravity. The usual Einstein-Hilbert 
action (with cosmological term) 



where the energy momentum tensor T^ v is obtained from the matter part 
of the Lagrangian, see §9.7 of [22] and §7.1.13 of [4T]. (Here we use the no- 
tation 70 for the (variable) cosmological constant, which is more frequently 
denoted by A or A in the cosmology and general relativity literature, but un- 
fortunately both of these letters are already assigned other meanings here.) 
In addition to these terms, where both G and 70 will be running with the 
energy scale A and depending on the free parameters /2 and fy, the asymp- 
totic expansion for the spectral action also delivers conformal gravity terms. 
Conformal gravity is considered an alternative to the usual form of general 
relativity, where the Einstein-Hilbert action is replaced by an action based 
on the Weyl curvature tensor 

Cx^uK = Rx^ivk — ^(gxuRfiK — 9fiuR\K + 9^kR\u) + -^{g\v9nK — 9Xk9^v)- 

This has the property of being conformally invariant, namely under a trans- 
formation of the form g^ u (x) 1— > f (x) 2 g ^{x) the Weyl tensor remains un- 
changed, Cx^i/k l— ► Cx^iuk- The action functional of conformal gravity is of 




minimally coupled to matter gives the Einstein field equations 

RV» - l g^R + lQ g^ = -8nGT^, 
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the form 

(4.1) a J C Xfl v K C x ^ K ^d 4 x, 

upon rewriting the above in terms of the Riemann curvature tensor 

and using the fact that R x ^kR X ^ k - AR^R^ + R 2 is a total divergence 
(see |37| ) one can rewrite the conformal action functional as 

2a J '{R^R"" - l -R 2 )^x 
which gives field equations 

1 

4a 

where 

W nu = 2C^ K - C^ Xuk R Xk = 

q9hu^ 2 R + -Vfj,V v R — -R{^R^u — -gfiuR) + R? p {-^9iivRfip — IRpupv)- 

In addition to the Weyl curvature tensor itself being conformally invariant, 
one can add to the conformal gravity action a coupling to a field (p. Under 
a conformal transformation gu V (x) i— > f(x) 2 g^ u (x), a field transforming like 
ip i— »• f~ 1( P gives 

(d^) 2 ' * r\{d^f + <^V>g/-2V^)VHog/), 

while the scalar curvature transforms like 

R ^ r 2 (R - (d - l)((d - 2) V M log /V log / + 2V^V M log /) 

so a non-minimal coupling of the field tp to gravity of the form 



2£o J R^ 2 ^gd A x 



is conformally invariant in dimension d = 4 if 2£o = 1/6. In addition to 
these terms a quartic potential 

Ao J <p A ^d 4 x 

also preserves conformal invariance. Thus, adding to the conformal action 

(4.1 ) terms of the form 

(4.2) i y {d^f^x - i y + A y 

maintains the conformal invariance. The conformal gravity action (4.1) 
with an additional non-minimal conformal coupling to another field ip in 
the form (4.2) is the basis of the Hoyle-Narlikar cosmologies, which were 
proposed as possible models for steady state cosmologies in [31 , with the 
field <p related to the Mach principle. In fact, the presence of the field 
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ip allows for a modification of the energy momentum tensor of the form 
T^ v i — ► T^ u - iw 2 {R» v - y^R) - g^Xoif*, which was used as "creation 
field" in steady state models. While the steady state cosmologies fail to 
account for major cosmological phenomena such as the background radia- 
tion, hence Hoyle-Narlikar cosmologies cannot be extrapolated towards the 
early universe, conformal gravity remains a valuable model (see [37] for a 
recent discussion). Within the noncommutative geometry model we will 
see below that one typically has a dominant Einstein-Hilbert action, and 
only at certain scales where the behavior of the running effective gravita- 
tional constant presents phase transitions one finds that the subdominant 
terms of conformal gravity become dominant. This gives rise to emergent 
Hoyle-Narlikar cosmologies, for which the problem of extrapolating towards 
earlier times does not arise, as they become suppressed by the dominant 
Einstein-Hilbert term away from the energy scale of the phase transition. 
In the noncommutative geometry model the role of the field non-minimally 
conformally coupled to gravity is played by the Higgs field. Variants of the 
model of [T7J, discussed for instance in [13], allow for the presence of a fur- 
ther scalar field a also conformally non-minimally coupled to R. Most of 
those arguments we describe in the following sections that are based on the 
coupling of the Higgs to gravity can be formulated also in terms of this other 
field a, though we will not explicitly mention it. 



4.2. Variable effective gravitational constant. According to the ex- 



pressions for the coefficients (2.12) of the terms in the asymptotic expansion 
of the spectral action, we see that this model has an Einstein-Hilbert term, 
where the usual coefficient 167rG J R^d^x, with G the Newton constant 
G ~ (10 19 GeV)~ 2 , with = 1.22086 x 10 19 GeV the Planck mass, is 

replaced by an effective gravitational constant of the form 

(4.3) G eff - K °- 3n 



8vr 192/ 2 A 2 - 2/„c(A) 
In this expresstion, the parameter /o is fixed by the unification condition 



(3.8) (we will use the first value for simplicity), the parameter / 2 is uncon- 
strained in and the function c(A) is determined by the renormalization 
group equations. 

Thus, one can best represent the effective gravitational constant of the 
model as a surface G e ff(A,/2), which is a function of the energy scale in 
the range A ew; ^ A ^ A-unif between the electroweak and the unification 
scales. It is in fact often preferable to consider the surface G~g(A,f2) that 
describes the inverse effective gravitational constant, since it is in this form 
that it appears in the coefficients of the action functional, and it does not 
have singularities. 

As an example of the form of the surface G~^(A, f'2), one sees in Figure [H] 
that, for sufficiently small /2, the surface exhibits a kink at the top see-saw 
scale. 
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As another example, if one requires in this model that the value of the 
effective gravitational constant at the lower end of the energy spectrum 
we are considering, that is, at the electroweak scale, already agrees with 
the usual Newton constant, this requires a very large fine tuning of the 
parameter /2, which is fixed to have the value f% = 7.31647 x 10 32 . Notice 
that this is just an example. There is no physical reason to assume in this 
model that at the time of the electroweak epoch of the early universe the 
effective gravitational constant would have to be already equal to that of 
the modern universe. In this example one then sees that the running of 
G~g becomes dominated entirely by the term 192/2 A 2 /(37r), with the term 
2/oc(A)/(37r) remaining all the while several orders of magnitude smaller. 
This is in contrast with the running of Figure [5] where fi is small, and the 
running of 2/oc(A)/(37r) comes to play a significant role, as one sees from 
the effect of the see-saw scales. 

4.3. Emergent Hoyle Narlikar cosmologies. We look more closely at 
cases with large f2- For simplicity we illustrate what happens in this range 
by focusing on the example we mentioned above where /2 is fixed so that 
G e ff(A ew ) = G, though qualitatively the results described in this section 
hold for a wider range of choices of /2 sufficiently large. We show that, in 
the case of spaces with R ~ 1, a phase transition happens at the top see- 
saw scale, where the dominant Einstein-Hilbert action is suppressed and the 
action of the system is dominated by a spontaneously arising Hoyle-Narlikar 
type cosmology. 

We first identify the dominant terms in the action, under the hypothesis 
of large fi- Then we identify the conditions under which these terms are 
suppressed and the remaining terms become dominant. 
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Proposition 4.1. For sufficiently large values of the parameter fa (for in- 
stance when G e g(A ew ) = G), and in the range of energies where the effective 
cosmological term is kept small by the choice of the parameter fa, the dom- 
inant terms in the expansion of the spectral action ( 2.11j ) are 



(4.4) A 2 ( JL j Ryfgd^x -ffij \H\ 2 ^-gd^ , 

for kq = Akq and flo = ixq/A. 

Proof. In this case, as we have seen already, the running of the effective grav- 
itational constant is dominated by the term G~g ~ 192/2A 2 /(37r). Similarly, 
in the quadratic term of the Higgs /x 2 ,, the term 2faA 2 / fa is dominant when 
fa is sufficiently large. In particular, in the example where G e g(A ew ) = G, 
the term 2faA 2 / fa satisfies 

2f2 f ew = 1.71627 x 10 36 , 
fa 

while the second term is several orders of magnitude smaller, 

-4^4 = -1-51201 x 10 27 , 
a{A ew ) 

even though the coefficient e varies more significantly than the other coeffi- 



cients of (2.8 ). 



We then proceed to estimate the remaining terms of (2.11 ) and show that 
they are all suppressed with respect to the dominant terms above, for this 
choice of fa. The parameters ao, To, £o are n °t running with A, and we can 
estimate them to be 

a 0.25916, r ~ 0.158376, 

while £0 = 1/12 remains fixed at the conformal coupling value. 

To estimate the running of the coefficient Ao, we propose an ansatz on 



how it is related to the running of A in the RGE (3.7) and to that of the co- 



efficients a and b. We also know that the boundary conditions at unification 

Ao|A=A u „ a/ = A(A uni /) 



of ( |2.12[ ) satisfy 

Tr 2 b(A unif ) 



faa 2 (A unif ) 



where in (3.7) one uses the boundary condition A(A un jj) = 1/2 as in [T]. 
Thus, we investigate here the possibility that the coefficient Ao runs with 
the RGE according to the relation 

(4.5) MA) = m 0£. 

With this ansatz for the running of Ao one finds that the value of Ao varies 
between \ (A ew ) = 0.229493 and A (A uni/ ) = 0.184711, along the curve 
shown in Figure k31 with a maximum at \a(A evJ ) and a local maximum at 
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the top see-saw scale with value 0.202167, and with a minimum near 2 x 10 7 
GeV of value 0.155271. We discuss this ansatz more in details in 94.5|below. 



Finally, for a given energy scale A it is always possible to eliminate the 
cosmological term at that energy scale by adjusting the coefficient (see 
{4.12). The vanishing condition for the effective cosmological constant at 



energy A is realized by the choice of 



(4.6) 



(4/ 2 A 2 c - f D) 
192A 4 



Thus, if at a given energy scale the coefficient f± is chosen so that the 
effective cosmological constant vanishes, then the only terms that remain as 
dominant terms in the action are those of (4.4). □ 



In this scenario then a new feature arises. Namely, the fact that one has 
comparable terms 



and 



1 

2kq 



96/ 2 A 2 
24vr 2 



2/ 2 A 2 



fo 



2.96525 x 10 32 A 2 



1.71627 x 10 32 A 2 



can lead to cancellations under suitable geometric hypotheses. 
Proposition 4.2. Consider a space with R ~ 1. Then, for values \H\ 



a/o/vr ; the term (4.4), which is the dominant term for fi sufficiently large, 



vanishes and is replaced by the sub-dominant 



(4.7) 



~\H\ 2 
a 



foe 
24vr 2 



R 



fo 



(e 



24' 



as the leading term in the formula (2.11) for the spectral action. Near the 



top see-saw scale, the term (4.7) has a discontinuity, where the dynamics 



becomes dominated by a Hoyle-Narlikar cosmology given by the remaining 



terms of (2.11 ). 
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Proof. As in [22], Corollary 1.219, we expand the Higgs field around \H\ 
Va/o/vr. Then, one can compare the two terms 

1 



(4.8) 



2kq 



R-ffi\H\ 



This identifies a value for the constant curvature 



R~ 

at which the dominant term 

96M 2 
24vr 2 



7T= 



R 



2/ 2 A 2 afo 
/( 



o 



of (4.8) vanishes, leaving the smaller terms to dominate the dynamics. One 



can estimate that this gives a value for the scalar curvature very close to 
one, R = 0.979907 ~ 1, if we use the value of a at unification energy and 



the first possible value of /o in (3.8) 



The next smaller term in (4.8) is then of the form 

e afo 



-\H\ 
a 



foe 
24vr 2 



R 



a 7r z 



foe 
24vr 2 



which gives (4.7). Near the top see-saw scale, at around 5.76405 x 10 14 GeV 
the term e — c/24 has a jump and a sign change due to the large jump of 
the coefficient e near the top see-saw scale (see Figure E| . At this phase 



transition what is left of the dynamics of (2.11) are the remaining terms 



One therefore sees an emergent behavior where near the phase transition of 
the top see-saw scale, of the following form. 

The coefficient f% can be chosen so that the cosmological term vanishes at 



this same top see-saw scale energy (see £4.12). The dynamics of the model 
is then dominated by the remaining terms, which recover a well known 
treatment of gauge and Higgs field in conformal gravity, as discussed for 
instance in §2.2 of [26 . According to these models, the conformally invariant 
action for the gauge and Higgs bosons is given by the terms 

1 



S c = «o / C tWfia C^^x + 



\DH\ z ^g-<Tx 



(4.9) 



-£o / R\H\ 2 ^d 4 x + \ J \H\ A ^d 4 



The additional topological term To f R* R* \fg d I x is non-dynamical and only 
contributes the Euler characteristic of the manifold. The action (4.9) is that 
of a Hoyle-Narlikar cosmology. □ 

In this model, the scalar curvature R, which is constant or near con- 
stant, provides an effective quadratic term for the Higgs and a correspond- 
ing symmetry breaking phenomenon, as observed in [26]. This produces 
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in turn a breaking of conformal symmetry, via the Higgs mechanism giv- 
ing mass to some of the gauge field, thus breaking conformal invariance. 
In this range, with a constant curvature R and in the absence of the qua- 
dratic term in the Higgs field is governed by a potential of the form 
Vr=i(H) = — £o\H\ 2 + Ao|-£f| 4 , which has a minimum at \H\ 2 = £o/(2\ ). 



We discuss in {4.7 below another instance of emergent Hoyle-Narlikar 
cosmology based on conformal gravity at a phase transition where the ef- 
fective gravitational constant undergoes a sign change and an antigravity 
regime arises. 

4.4. Effects on gravitational waves. The fact that we have a variable 
gravitational constant in the model has detectable effects on phenomena like 
gravitational waves whose propagation depends on the value of the gravita- 
tional constant. 

Under the assumptions that the remaining terms in the asymptotic ex- 
pansion of the spectral action are negligible with respect to the dominant 



(4.4), and further that \H\ ~ 0, so that only the Einstein-Hilbert term dom- 



inates, one can show as in [38 that the equations of motion for (2.11 ) reduce 
to just 

In [38] the authors conclude from this that, in the isotropic case, the non- 
commutative geometry model has no effect on the gravitational waves that 
distinguish it from the usual Einstein-Hilbert cosmology. However, in fact, 
even in this case one can find detectable effects on the gravitational waves 
that distinguish the noncommutative geometry model from the ordinary 
case of general relativity, because of the running of the effective gravita- 
tional constant. 

We consider here two different scenario for the time-energy relation, one 
which will be relevant close to the electroweak scale ([30]), where A ~ i -1 / 2 , 
and the other that refers to the inflationary period, with A ~ e~ at . We 
show that in both cases the behavior of the gravitational waves differs from 
the behavior, with the same time-energy conversion, of the solutions in the 
classical case, thus detecting the presence of noncommutative geometry. 

Proposition 4.3. In between the unification and the electroweak scale the 
gravitational waves propagate according to the A- dependent equation 

a\ 2 If. I a\ ■ ~\ 12vr 2 



< 410 > - 3 {a) + -2{ i {a) h+2h ) = ^-MA) T '"- 

Upon rewriting the energy variable A as a function of time through A 
l/a(t), one obtains 



(4.11) 
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in the radiation dominated era where A ~ i" 1 / 2 , and 



2 r , 



(4.12) 



h + 2ah — 3a 



12vr 2 r 00 
96f 2 e~ 2at - Me 1 "*) ' 

in the inflationary epoch where A ~ e~ at . 

Proof. For a metric of the form 



(4.13) 



9 to, = a{tf 



-1 







5ij + hij(x) 



one separates the perturbation hij into a trace and traceless part, and the 
gravitational waves are then governed by the Friedmann equation, which 
gives 

2 



(4.14) 



+ 



1 



h + 2h 



This equation is formally the same as the usual equation for the gravita- 
tional waves, up to replacing acq for 87rG, as remarked in [38]. However, the 
dependence of Kg on the energy scale A leads to the result (4.10). 

The change of variable between energy and time, for a cosmology of the 
form (4.13) is given by A = l/a(t). Thus, we can write the equation (4.14) 

side 

a(t)' 



by expressing the right hand side also as a function of time in the form 

2 



(4.15) 



a(t) 
a(t) 



+ 



1 



o(t) 



h(t) + 2h(t) 



96/ 2 



127r 2 r 00 

J Mwj- 



(a(t)) 



In the radiation dominated era, where the function a(t) behaves like a(t) 
t 1 / 2 we find 

1 



a 
a 



-t 



-i 



which gives 



-3 



+ 



h + 2h 



3 

h + t^h- -t~ 2 . 
4 



This gives (4.11). In the inflationary era, where the function a(t) behaves 



exponentially a(t) = e , one obtains instead a/a = a and 

2 



+ 



1 



h + 2h) = h + 2ah-3a. 



which gives (4.12). 



□ 



For a choice of the parameter fi sufficiently large (such as the one that 
gives G e g(A ew ) = G) for which Kq(A) ~ Kq/A 2 we then obtain the following 
explicit solutions to the equation (4.10). 
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Proposition 4.4. We consider the case with the parameter f 2 sufficiently 
large (for instance in the example where G e s(A ew ) = G). In the radiation 
dominated era, where the metric (4.13) has a(t) ~ t x l 2 , and the energy-time 
relation is given by A = 1 /a(t) the equation (|4.11j) has solutions of the form 



(4.16) 



hit) 



^t* + B + A\o g {t) + \^ g {t) 2 



In the inflationary epoch, for a metric of the form (4.13) with a(t) 
for some a > 0, the equation (4.12) has solutions of the form 

(4.17) 



~ e 



at 



h(x) 



3tt 2 T ( 



I92f 2 a 2 



00 e 2at 



3a A n.. „ 
+ ^t+—e~ 2at + B. 
2 2a 



Proof. We assume that the parameter f 2 is sufficiently large. For example, 
we take the case where f 2 is fine tuned so as to have an agreement between 
the value at the electroweak scale of the effective gravitational constant and 
the usual Newton constant, Kq(A £W ) = 87rG. One then finds at high energies 
a different effective Newton constant, which behaves like 

2 12 ^A-2 

96/ 2 

which then gives a modified gravitational waves equation 



(4.18) 



,-.2 



A 2 



Too, 



with h 2 = 



We look first at the radiation dominated case where a{t) = i 1//2 . In the 
interval of energies that we are considering, this has relevance close to the 



electroweak scale, see |30J . We then have equation (|4.11|) in the form 
(4.19) 



12vr 2 Tf 



00 



Assuming Too constant, the general solution of 



t~ r h + t~ 2 h- -t~ 



C 



is of the form 



h(t) = ^t 3 + B + Alog(t) + ^log(t) 2 , 



for arbitrary integration constants A and B. With C = (127r 2 Too)/(96/2), 
this gives (4.16). 

In the inflationary epoch where one has a{t) = e at , for some a > 0, one 
can write the equation (4.12) in the form 



h + 2ah - 3a 2 = e 2at 



127r 2 r 00 
96/ 2 



EARLY UNIVERSE MODELS FROM NONCOMMUTATIVE GEOMETRY 



27 



Again assuming Too constant, the general solution of 

h + 2ah - 3q 2 = Ce 2at 
with a constant C is of the form 

h(+\ — _i_ Jl_ c 2crf A -2at , R 

for arbitrary integration constants A and B. With C = (12-7r 2 Too)/(96/2) as 



above, this gives (4.17). □ 



The behavior of the solutions (4.16) and (4.17) should be compared with 



the analogous equations in the ordinary case where of the equation (4.14) 



with Kq = 87rG Too independent of the energy scale A and equal to the ordi- 
nary Newton constant. In this case, one obtains, in the radiation dominated 
case with a(t) = t 1 / 2 the equation 

3 

h + t' l h- -r 2 = 8vrGToo, 
which has general solution 

h(t) = 2itGT 00 t 2 + B + Alog(t) + I log(i) 2 , 



which differs from (4.16) for the presence of a quadratic instead of cubic 
term. Similarly, in the case of the inflationary epoch, where one has a{t) = 
e at , in the ordinary case one has the equation 



h + 2ah — 3a 2 = 8^(7 Too, 



which has general solution 



( 47rGT00 + ^)t + ^e-^ + B, 
a 2 2a 



which differs from what we have in (4.17) by the presence of an additional 



linear term instead of an exponential term. 

Similar examples with different forms of the factor a(t) in the metric 



(4.13) can easily be derived in the same way. Choices of the parameter 
/2 for which the term c(A) cannot be neglected will give rise to varying 
behaviors of the equations both in the radiation dominated era and during 
inflation. However, in those cases the equations cannot be integrated exactly 
so we cannot exhibit explicit solutions. 

4.5. The Ao-ansatz and the Higgs. We mention here briefly another con- 



sequence of the ansatz (4.5) on the running of the coefficient Ao(A) between 
the electroweak and the unification scales. This running is different from 
the one used in [T7] to derive the Higgs mass estimate. In fact, the RGE 
equations themselves are different, since in [T7] one only considers the RGE 
for the minimal standard model and the boundary conditions at unifications 
are also significantly different from the ones used in [I] that we use here (see 
the second appendix for a more detailed discussion of the boundary condi- 
tions in [T7| and [1]). Also in [T7j the coefficient Aq is assumed to run as 
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the A in the minimal standard model, with only the boundary condition at 
unification relating it to the values of a and b (which are in turn different 
from the values at unification according to the RGEs of [1 ). Thus, one can 



check how adopting the ansatz (4.5) for Ao(A) together with the boundary 
conditions of [1 for A, a and b affects the estimate for the Higgs mass. In 
|17] one obtains a heavy Higgs at around 170 GeV, by the value 

V2A— ~ 170 GeV, 
9 

where A is the low energy limit of the RGE flow in the minimal Standard 
Model for the coefficient Ao and where 2M/g ~ 246 GeV is the Higgs vac- 
uum. If we replace the running used in [17 with the running of Aq(A) of 



(4.5) with the boundary condition of [T], the same estimate would deliver a 



much lower value 



g y f a 2 (A ew ) g 

This looks potentially interesting in view of the fact that the projected 
window of exclusion for the Higgs mass in [16] starts at 158 GeV (see also 
|32] ) . This gives only a first possible indication that a more detailed analysis 
of the RGEs for the standard model with Majorana mass terms, as in [1] , 
and a careful discussion of the boundary conditions at unification and of 
the running of the coefficients in the asymptotic expansion of the spectral 
action may yield a wider spectrum of possible behaviors for the Higgs field 
within these noncommutative geometry models. This topic deserves more 
careful consideration that is beyond the main focus of the present paper. 

4.6. Antigravity in the early universe. Cosmological model exhibiting 
a sign change in the effective Newton constant in the early universe, due to 
the interactions of gravity and matter, were studied for instance in [35] , [H] 
or [IE]. Those models of antigravity are based on the presence of a non- 
minimal conformal coupling of gravity to another field, with a Lagrangian 
of the form 

C = I 16^G^ + 1^ / + 4^), 

where the last term contains the kinetic and potential terms for the field cp 
and its interactions with other fields A, ip. The conformal coupling of R and 
<p gives rise to an effective gravitational constant of the form 

G-fj 1 = G- 1 - A - V- 

where tp is treated as a constant, which is estimated in f32] in terms of the ex- 
cess of neutrino over antineutrino density, the quartic interaction coefficient 
A of p and the Weinberg angle. A decrease in G e g produces a correspond- 
ing increase in the Planck density (see, however, the criticism to this model 
discussed in [H]). Antigravity sectors with negative effective gravitational 
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Figure 7. An example of transition to negative gravity in 
the running of G~$(A) with G cS (A unif , f 2 ) = G. 

constant as in (35] were recently considered within various approaches based 
on extra dimensions and brane world models, see for instance [29], [28], or 
from the point of view of moduli in heterotic superstring theory, as in [12] . 

Within the noncommutative geometry model it is also possible to find sce- 
narios where one has antigravity in the early universe. The same mechanism 
proposed in [35] can be reproduced within the NCG model, due to the pres- 
ence of the conformal coupling of the Higgs field to gravity. However, there 
is another, independent mechanism that can also produce a sign change of 
the Newton constant between the unification and electroweak cosmological 
phase transitions and which is only due to the running of the effective New- 
ton constant with the RGE equations of the particle physics content of the 
model. In fact, there are choices of the parameter f 2 of the model for which 
the effective gravitational constant undergoes a sign change. 

An example of this behavior is obtained if one chooses the value of the 
effective gravitational constant G c s to be equal to the Newton constant 
at unification scale, G e s(A un if, f 2 ) = G, and then runs it down with the 
RGE equations. Notice that the assumption G e s(A un if, f 2 ) = G is the same 
that was proposed in [T7] , but due to the different RGE analysis considered 
here, the behavior we describe now is different from the one projected in [T7] . 
With the estimate used in Lemma 5.2 of [17] for the Majorana mass terms at 
unification (under the assumptions of flat space and negligible Higgs vacuum 
expectation) setting the effective gravitational constant of the model equal 
to the Newton constant at unification energy only requires f 2 of the order 
of at most f 2 ~ 10 2 , while using the boundary conditions of [1] one find a 
larger value of f 2 . In fact, we see that setting 

(4.20) f 2 ~ 18291.3 

gives G eS (A unif ) = G = (1.22086 x lO 19 )" 2 . 

One sees then, as in Figure [7]that the resulting G e g(A, f 2 ), for this choice 
of f 2 has a sign change at around 1.3 x 10 12 GeV. Thus, with these boundary 
conditions one finds an example of a regime of negative gravity in the early 
universe. Other possible choices of f 2 lead to similar examples. 
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4.7. The conformal gravity regime. In this scenario, it is especially in- 
teresting to see what happens near the energy scale of 1.3 x 10 12 GeV, where 
G~g vanishes, as in Figure [TJ This gives another example of an emergent 
conformal gravity regime at a phase transition (here the change from positive 
to negative gravity) of the system. 

Proposition 4.5. Let the parameter fi be chosen so that the inverse ef- 
fective gravitational constant G~^(A) has a zero at some A = Ao- Assume 
the vanishing of the topological term and suppose that the parameter f$ is 
chosen so that the effective cosmological constant also vanishes at Aq. Then 



near Ao and for \H\ 2 sufficiently small the dynamics of (2.11) is dominated 
by pure conformal gravity. 

Proof. At the singularity A s i ng for G e s, assuming the vanishing of the topo- 
logical term, the terms that remain in the bosonic part of the action are 
the cosmological term, the conformal gravity term with the Weyl curvature 
tensor, and the Higgs and gauge bosons terms. If the Higgs field is suffi- 
ciently near the H = vacuum, and the parameter f± is chosen so that the 
cosmological term also vanishes at the same scale A s i ng , one finds that what 
remains of the bosonic action is just the conformal gravity action 

(4.21) S(A smg ) = a J C, lvpu C^Jgd A x 

with ao = — 3/o/(107r 2 ) ~ —0.25916, and with the additional weakly coupled 
term of the gauge bosons. Thus, in this scenario, when running down the 
coefficients of the bosonic spectral action from unification scale towards the 
electroweak scale, a singularity of G e g occurs at an intermediate scale of 
1.3 x 10 12 GeV. At the singularity, if the cosmological term also vanishes and 
the Higgs contribution is sufficiently small, the model becomes dominated 



by a conformal gravity action (4.21). 



In this regime the equations of motion will then be of the form 
(4.22) 2C^ UK - C" Xuk R Xk = — 

,A„K 4 a[) 

as in (188) of [37J. Thus one has in this regime pure conformal gravity with 
an effective gravitational constant of see §8.7 of |37j. □ 

The emergence of a conformal gravity regime was also observed in [38 , 
though not in terms of renormalization group analysis, but near the special 
value of the Higgs field \H\ — > \/6/ko- 

4.8. Gravity balls. We now analyze more carefully the second mechanism 
that produces negative gravity besides the running of G e ff(A), namely the 
non-minimal conformal coupling to the Higgs field. This will reproduce in 
this model a scenario similar to that of |35| . 

Let us assume for simplicity that the parameter fi is chosen so that 
GW(A, f-i) > for all A ew < A < A uni j. We show that, even in this case, 
it is possible to have regions of negative gravity, due to the coupling to the 
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Higgs. These behave like the gravity balls and non-topological solitons of 
[36] , jl3] , but with a more elaborate behavior coming from the fact that the 
underlying gravitational constant is also changing with A according to the 
RGE flow. 

Proposition 4.6. Let f 2 be assigned so that the effective gravitational con- 
stant satisfies G e g(A, f 2 ) > for allA ew < A < A uni f. Then negative gravity 
regions, with \H\ near \H\ 2 ~ /Xq/(2Ao), arise in the range of energies A such 
that 



(4.23) 
for 

(4.24) 

and 

(4.25) 



e H (Aj 2 )>e G (Aj 2 ), 



MA, /a) 



e G (Aj 2 ) 



(2/ 2 A 2 q(A) - / e(A))a(A) 
7r 2 A(A)b(A) 

192/ 2 A 2 - 2/ c(A) 
4tt 2 



Proof. The presence of the conformal coupling term 

-1 J R\H\ 2 ^gd 4 x 



in the normalized asymptotic formula for the spectral action (2.11) means 
that, in regions with nearly constant \H\ 2 , the effective gravitational con- 
stant of the model is further modified to give 



(4.26) 



1 



f Geffl^l 2 ' 



This is the same mechanism used in [35J for negative gravity models. This 



means that, assuming that G e g(A) > for all A e 
have 

' G cS ,h < for \H\ 2 > 
G eSH >0 for \H\ 2 < 



< A < A un if, one will 



4vrG cff (A) ' 
3 



47rGW(A) 

This means, for instance, that in the presence of an unstable equilibrium at 
|iT| = and a stable equilibrium at \H\ 2 = v 2 satisfying v 2 > ^g^A) > one 
can have gravity balls near zeros of the field \H\ 2 , where gravity behaves in 
the usual attractive way, inside a larger scale negative gravity corresponding 
to the true equilibrium \H\ 2 = v 2 . 



In the action (2.11 ), the Higgs field has a quartic potential given by 



-Mo 



j \H\ 2 ^d 4 x + \ Q j \H\ 4 ^d 4 x, 



:!2 



MATILDE MARCOLLI AND ELENA PIERPAOLI 




1x10' 



Figure 8. The surface ^f/(A,/ 2 ) in the range 10 14 < A < 
10 15 GeV, around the top see-saw scale, for 10~ 4 < / 2 < 1. 

which has a minimum at /iQ/(2Ao). Thus, to identify the negative gravity 



regime we need to compare the running of (4.24), 

(2/ 2 A 2 a(A) - / e(A))o(A) 



^(A, f 2 ) 



2A 



(A) 



/2 A 2 
1 fo 

A(A) 



e(A) 
a(A) 



/oa 2 (A) 



7T 2 A(A)b(A) 



where we again used the ansatz (4.5) on the running of the coefficient Ao(A), 

on ([425}, 

3 192/ 2 A 2 - 2/ c(A) 192/ 2 A 2 - 2/ c(A) 



and the running of the function (|4.25|), 
3 



fc(A,/ 2 ) := 



4tt 



3tt 



4^ 2 



4vrG cff (A) 

Thus, for different possible values of the parameter / 2 , a negative gravity 
regime G e s h < and gravitational balls are possible in the range where 
^(A,/ 2 ) >£ G {A,f 2 ). □ 

An example of such a transition to a negative gravity region is illustrated 
in Figure |9] where the surface ^#(A, / 2 ) near the top see-saw scale behaves 
as in Figure [8j 

4.9. Primordial black holes. The possibility of primordial black holes 
(PBHs) in the early universe was originally suggested by Zeldovich and 
Novikov in the late 1960s (see [SS] for a survey). They originate from the 
collapse of overdense regions, as well as from other mechanisms such as phase 
transitions in the early universe, cosmic loops and strings, or inflationary 
reheating. The consequences on primordial black holes of the running of the 
gravitational constant in the early universe have been analyzed, for example, 
in In models where the gravitational constant in the early universe 
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Figure 9. Transition to a negative gravity regime where 
M A >/2) > ^g(A,/ 2 ) in the region 10 12 < A < 10 15 with 

/2 = 1- 



may be different from the value it has in the modern universe, PBHs, whose 
existence is conjectured but has not been presently confirmed, are seen as 
a possible source of information about the changing gravitational constant. 
In fact, the mass loss rate due to evaporation via Hawking radiation and 
the Hawking temperature of primordial black holes depends on the value of 
the Newton constant so that the evolution of such black holes depends on 
the change in the Newton constant. The evaporation of primordial black 
holes is often proposed as a mechanism underlying 7-ray bursts (see for 
instance [6], [9]). An especially interesting question regarging PBHs is that 
of gravitational memory, as described in [5]. It is usually assumed that two 
possible scenarios for the evolution of primordial black holes can happen: 
one where the evolution follows the changing gravitational constant and one 
with the possibility of "gravitational memory" , namely where the evolution 
of the PBH is determined by an effective gravitational constant different 
from the one of surrounding space. This latter phenomenon can arise if the 
model has the possibility of having regions, as in the case of the gravity 
balls discussed in §4.8| above, where the effective gravitational constant has 
a value different from the one of surrounding space. 

We discuss here the effect on primordial black holes of the running of 
the effective gravitational constant in the noncommutative geometry model. 
For primordial black holes that formed at a time in the early universe when 
the gravitational constant G e g (i) was different from the one of the modern 
universe and whose evolution in time reflects the corresponding evolution of 
the gravitational constant, the black hole "adjusts its size" to the changing 
G e g(t) according to the equation 



(4.27) 
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where Ai is the mass of the primordial black hole. Correspondingly the 
temperature varies with the changing gravitational constant as 

-i 



(4.28) 



(8irG cS (t)M(t)Y 



In the case with gravitational memory, one can have a black hole that evolves 
according to a gravitational constant that is different to the one of the 
surrounding space. In our setting, given the scenario described in £4.8, this 
second case occurs with the equation (4.27) replaced by 

dM{t) 



(4.29) 



dt 



-(G eS , H (t)M(t)Y 



with Geg^H as in (4.26). 



Proposition 4.7. In the radiation dominated era, for a cosmology with 
metric tensor of the form (4.13), the evaporation of primordial black holes 
by Hawking radiation in the NCG model is given by 

(4.30) M(A,f 2 )-- 



! M 3 (A l 



2 
3^2 



(192/ 2 x 2 - 2/ c(x)) 2 



dx, 



in the case of PBHs without gravitational memory, while PBHs with gravi- 
tational memory evaporate according to 



(4.31) 

where 



M(A,f 2 ) = ^jM 3 (Ai 

192/ 2 A 2 - 2/ c(A) 
3vr 



2 

37T 2 



in (1 



fG eff (*)|g| 2 ) 2 
r3n Zf^a ax > 



x 3 G eS (x) 



[I 



4fG eff (A)|#| 2 ) 



G cS (A) 
4 (2/ 2 A 2 a(A) 



/oe(A))a(A) 



3?r 



A(A)b(A) 



for \H\ 2 ~ /u|/(2A Q ). Here A i n is the energy scale at which the radiation 
dominated phase begins, that is, where one starts to have a(t) = i 1 / 2 . 



Proof. For a metric of the form (4.13), in the radiation dominated era one 
has a(t) ~ i 1//2 : , hence the energy-time change of variables A = \/a(t) is 
of the form A = t -1 / 2 . Thus, the equation (4.27) can be rewritten in the 
variable A in the form 

dM(A) 2 



(4 " 32) dA A3(G eff (A)^(A)) 2 ' 

since dt/dA = -2A~ 3 . We look at the case of a primordial black hole that 
has a given mass «M(Aj n ) at the initial Aj n , and we look at the evolution of 
«M(A) between Aj n and electroweak scale. The equation (4.32), written in 
the form 



M{A) 2 dM{A) 



dA 
! A 3 G 2 ff (A) 
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gives (4.30). The case with gravitational memory is obtained similarly by 
replacing the effective gravitational constant G e g(A) with the one locally 
modified by the interaction with the Higgs field, G e ff #(A). □ 

The expression for the evaporation by Hawking radiation simplifies in the 
cases where is sufficiently large that the term 192/2 A 2 dominates over 
2/oc(A). This is the case, for example, when G e g(A ew ) = G and is chosen 
accordingly. 

Corollary 4.8. In the case where /2 is sufficiently large that the term 
2/oc(A) is negligible with respect to I92/2A 2 for all A ew < A < A,, n; the 
evaporation law for a primordial black hole that forms at Aj n gives a bound 
on its mass by (Aj n 64v / 6/2/ 7r ) 2 ^ 3 - 

Proof. Under the assumption that I92/2A 2 dominates over 2/oc(A), the right 
hand side of (4.32) can be approximated by the dominant term of G e g(A) 
which gives 

(4.33) M 2 dM=2^^j AdA. 



Let Ao denote the scale at which the PBH evaporates. Then (4.33) gives 
^(A) = 2(^M) 2 (A 2 -A 2 ). 

This sets a bound to the mass of a primordial black hole formed at A = A». 
from the condition that 

Ao = A 2 „ — ( — — — -=- \ Mu ni f 



which gives 



Mir, < 



(64/ 2 V6) 



A m 64/ 2 V6 V /3 



□ 



We can analyze similarly the equation of the Hawking radiation for PBHs 
during the inflationary epoch when one has a(t) = e at for some a > 0. We 
obtain the following result. 

Proposition 4.9. In the inflationary epoch, for a cosmology with metric 
tensor of the form (4.13), the evaporation of primordial black holes by Hawk- 
ing radiation in the NCG model is given by 
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in the case of PBHs without gravitational memory, while PBHs with gravi- 
tational memory evaporate according to 

<« 5 > ma. h) - f^> - ^ f; (1 - !°;tT' 2)2 

where Aj n is the energy scale at which the inflationary behavior a(t) = e at 
begins. 

Proof. The argument is completely analogous to the previous case. Here the 
energy-time change of variables A = l/a(t) gives t = —a -1 log A. Thus, in 
the energy variable the equation ( 4.27| ) becomes 

aA(G eff (A)) 2 

for the case without gravitational memory, or the same equation with G e s (A) 
replaced by G e fr ) H(A) in the case with gravitational memory. This gives 



(4.34) and (4.35). □ 



4.10. Higgs based slow- roll inflation. Recently, a mechanism for infla- 
tion within the minimal standard model physics was proposed in [24 . It is 
based on the presence of a non-minimal coupling of the Higgs field to gravity 
of the form 



-Co / R\H\ 2 ^d 4 x 



as we have in the asymptotic expansion of the spectral action in the non- 
commutative geometry model, but where the value of Co is n °t set equal to 
the conformal coupling Co = 1/12, but is subject to running with the RGE 
flow. In [21] the running of Co is governed by the beta function given in |8 , 
which in our notation we can write approximately as 

(4.36) 16^/%, (A) = (-12Co(A) + 1) F(Y U , Y d , Y„, Y e , M, g u g 2 , g 3 ,X), 

where the function F of the running parameters of the model is computed 
explicitly in [8]. In [21] this running is only considered within the minimal 
standard model, without the right handed neutrinos and Majorana mass 
terms M, but their argument can be adapted to this extension of the stan- 
dard model, since the general derivation of the running of Co in [BJ applies in 
greater generality. In this variable Co scenario, the dimensionless quantity 
that governs inflation is tp = \/£o\H\/mp, where mp is the reduced Planck 
mass, which in our notation is mp = 1/kq. The inflationary period corre- 
sponds in [21] to the large values tp >> 1, the end of the inflation to the 
values ip ~ 1 and the low energy regime to ip « 1. 

At present we do not know whether a modification of the noncommutative 
geometry model that allows for a variable Co (A), different from the conformal 
coupling Co = 1/12 is possible within the costraints of the model, but we 
show here that, even in the case where Co is constant in A and equal to the 
conformal Co = 1/12, the noncommutative geometry model still allows for 
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Figure 10. Plot of V E (H)/V E against i/j. 

a similar inflation mechanism to occur, through the running of the effective 
gravitational constant. 

Proposition 4.10. A Higgs based slow roll inflation scenario arises in the 
NCG model with parameter 

V>(A) 2 = £o^(A)|F| 2 

and with potential 

(4.37) V E (H) A(,|//|1 



(i + £o«g|# 



2\2 ' 



Proof. In the noncommutative geometry model, the coefficient kq(A) is run- 
ning with A according to 



«o(A) 



12vr 2 



96/ 2 A2 - / c(A) • 
Then the parameter that controls inflation is given by 

If one then proceeds as in [24J, one finds that In the Einstein metric 
9% = f(H)gfj, v , for f(H) = 1 + £oko\H\ 2 the Higgs potential becomes of 



the form ( 4.37| ). In the range where ip » 1 this approaches the constant 



function (constant in H but not in A) 

A (A) A(A)b(A)(96/ 2 A 2 -/ c(A)) 2 



V E 



4£ 2 (A)^(A) " 4/ a 2 (A) 



where we used again the ansatz (4.5 ) for the running of the coefficient Ao(A), 
while at low values ip « 1 the potential is well approximated by the usual 
quartic potential Ve(H) ~ Ao|-£f| 4 . □ 

The asymptotic value Ve in turn depends on the energy scale A and 
different behaviors are possible upon changing the values of the parameter 



fi of the model as the examples in Figure 11 illustrate. One can see the 
effect of the top see saw-scale on the running. 
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4xl0 14 6xl0 14 8xl0 14 1x10" 



Figure 1 1 . The behavior of the asymptotic value Ve of the 
potential Ve(H), plotted above for f 2 = 10~ 4 , / 2 = 10~ 3 , 
and h = 1/10 and for 10 14 < A < 10 15 



Notice that a potential of the form as in Figure 10 which is especially 
suitable for slow roll inflationary models, also arises naturally in the NCG 
setting from non-perturbative effects in the spectral action, as in [16 . We 
plan to return in future work to describe how the non-perturbative approach 
can be applied to cosmological models beyond the very early universe. 

4.11. Spectral index and tensor to scalar ratio. In a slow roll inflation 
model, the first and second derivatives of the potential Ve together with a 
change of variable in the field that brings the action into a canonical form, 
determine the first and second slow-roll parameters e and r), see [Mj §3.4 
and also the analysis in [24 . We analyze here the form of the slow-roll 
parameters in the noncommutative geometry model. 



We write the potential Ve of (4.37) in the form 

A s 4 



(4.38) V E (s) 



(l + £ ^ 2 ) 2 ' 



We will consider derivatives of this potential as a function of the variable s. 
We also introduce the expression 

(4.39) C(.)-= „„ ' + 3 (2& *' 2 



which corresponds to the (da / dcj)) 2 of equation (4) of [23] that gives the 
change of variable in the field that runs the inflation that puts the action in 
a canonical form. 
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The two slow-roll parameters are then given by the expressions 

2k 2 \ V E {s) 



which correspond to the equations (6) and (7) of [23]. We then have the 
following result. 



Proposition 4.11. The first slow roll coefficient is of the form 

(442) £(S) = S 2 + &^ 1 (l + («g) 2 )- 4 

while the second slow-roll coefficient is given by 

, 4 43) (s) = 8(3 + £o^ 2 (l - 2^4{s 2 + 12«g(-l + fr)/^ 2 )))) 

K 2( S + ^(1 + (^) 2 )S 3 ) 2 

Proof. The derivatives of the potential are 

, 4AoCoKoS 5 4A s 3 

V E(s) = - TTT1- -2^3 + 



(1+^ 2 ) 3 (1+^ 2 S 2 ) 2 

t r// / v _ 24A £o(*o) 2 s 6 36A £ ^ 4 , 12A o* 2 



(1 + eo^s 2 ) 4 (i + ^ 2 ) 3 (i + ^ 2 ) 2 ' 

One also has 

C(s )-3/2 = >/8(l+ft«fc 2 ) 3 



V(i+^(i + i2^)s 2 ; 



3 



A(7( s )i/2 = 6«§s(l + go^(s 2 + 12k 2 (-1 + go^s 2 ))) 



^ (l + £o^ 2 ) 2 ^2(1 + eo«§(l + 12£o(^) 2 )* 2 

One then computes e(s) and 77(5) directly from (4.40) and ( |4.41[ ). □ 

Notice that both coefficients depend on the energy scale A and on the 
parameters fi and /o through the single parameter 

«o(A) 



96/ 2 A 2 - / c(A) ' 

while no dependence through Ao(A) remains and in our model £0 = 1/12. 

The slow-roll coefficients provide expressions for the spectral index n s and 
the tensor to scalar ratio r, which can be directly compared with cosmologi- 
cal data. Thus, a more detailed analysis of the behavior of these as functions 
of A and of the parameter fi of the model will give an exclusion curve for 
the parameter. We will provide a more detailed analysis elsewhere, but for 
the purpose of the present paper we derive the corresponding expression one 
has in this model for the spectral index and the tensor to scalar ratio. 
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Proposition 4.12. The spectral index is of the form 
(AAA) 

32(216 + n 2 {6s 2 - «g(432 + 12«g(2 + 3(/tg) 2 )s 2 + (1 + (/^) 2 )s 4 ))) 

while the tensor to scalar ratio if given by 

(4.45) r = 2 y . 

, 2 + f§(l + (K 2 )2) S 4 

Proof. One knows that the spectral index and the tensor to scalar ratio are 
related to the slow-roll coefficients by (see [M] §7.5.2 and §7.6 and [24] ) 

(4.46) n s = 1 - 6e + 2r?, and r = 16e. 



The result then follows directly from Proposition 4.11. □ 



Once again, these parameters have a dependence on the choice of the free 
parameter fi of the model and a scaling behavior with the energy A, which 
depends on the running of Kq(A). Since the spectral index and the tensor to 
scalar ratio are heavily constrained by cosmological data from the WMAP 
combined with baryon acoustic oscillations and supernovae data, this pro- 
vides a way to impose realistic constraints on the parameter fi of the model 
based on direct confrontation with data of cosmological observations. We 
will provide a detailed analysis of the constraints imposed on the noncom- 
mutative geometry model by the spectral index and tensor to scalar ratio in 
a separate paper. Since the parameters n s and r are also providing infor- 
mation on the gravitational waves (see [33] §7.7) a more detailed analysis of 
their behavior and dependence on the parameters of the model will give us 
a better understanding of the effects on gravitational waves of the presence 
of noncommutativity and may explain in a different way the amplification 
phenomena in the propagation of gravitational waves that we observed in 



4.12. Variable effective cosmological constant. The relation between 
particle physics and the cosmological constant, through the contribution of 
the quantum vacua of fields, is well known since the seminal work of Zel- 
dovich [19] . The cosmological constant problem is the question of reconciling 
a very large value predicted by particle physics with a near zero value that 
conforms to the observations of cosmology. Among the proposed solutions 
to this problem are various models, starting with [3£J], with a varying ef- 
fective cosmological constant, which would allow for a large cosmological 
constant in the very early universe, whose effect of negative pressure can 
overcome the attractive nature of gravity and result in accelerated expan- 
sion, and then a decay of the cosmological constant to zero (see also [ID] for 
a more recent treatment of variable cosmological constant models). Often 
the effective cosmological constant is produced via a non-minimal coupling 
of gravity to another field, as in [25] , similarly to what one does in the case 
of an effective gravitational constant. 
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In the present model, one can recover the same mechanism of [25] via the 
non-minimal coupling to the Higgs field, but additionally one has a running 
of the effective cosmological constant 70(A) which already by itself may 
produce the desired effect of decaying cosmological constant. We illustrate 
in this section an example of how different choices of the parameter fa, 
for a fixed choice of fa, generate different possible decay behaviors of the 
cosmological constant. These can then be combined with the effect produced 
by the non-minimal coupling with the Higgs field, which behaves differently 
here than in the case originally analyzed by [25]. This still does not resolve 
the fine tuning problem, of course, because we are trading the fine tuning 
of the cosmological constant for the tuning of the parameters fa and fa of 
the model, but the fact that these parameters have a geometric meaning in 
terms of the spectral action functional may suggest geometric constraints. 



As we have seen in the proof of Proposition 4.1 above, one can impose the 



vanishing of the effective cosmological constant 70 at a given energy scale 



A by fixing the parameter fa equal to the value given in (4.6). This means 
that, when the effective gravitational constant of the model varies as in the 
effective gravitational constant surface G c s(A, fa), depending on the value 
assigned to the parameter fa, there is an associated surface that determines 
the value of the parameter fa that gives a vanishing effective cosmological 
constant. This is defined by the equation 

(4.47) tt 2 7o(A, fa, fa) = 48/ 4 A 4 - / 2 A 2 c(A) + i/ D(A) = 0. 

The solutions to this equation determine the surface 

/ 2 A 2 c(A) - i/ 0(A) 



(4.48) fa(A,fa) 



48A 4 



To illustrate the different possible behaviors of the system when impos- 
ing the vanishing of the effective cosmological constant at different possible 
energy scales A, we look at the examples where one imposes a vanishing 
condition at one of the two ends of the interval of energies considered, that 
is, for A = A ew or for A = A un if. For simplicity we also give here an explicit 
example for a fixed assigned value of fa, which first fixes the underlying 
G e ff(A). we choose, as above, the particular example where fa is chosen 
to satisfy G e g(A ew ) = G. Under these conditions we impose the vanishing 
of the effective cosmological constant at either A ew or A un if through the 
constraint 

(4.49) 7o{A ew , fa, fa) = or 7o(A Mni/ , fa, fa) = 

which fixes the value of fa = fa(A ew , fa). Then the behavior at higher 
energies of the effective cosmological constant is given by the surface 

(4.50) lo{A,fa,fa{A ew ,fa)) or 7o (A, fa, fa{A unif , fa)) 
as a function of the energy A and the parameter fa. 
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Figure 12. Running of 70(A) under the assumptions of van- 
ishing 7o(A ew ) = or 7o(A mi/ ) = 0. 



We find two very different behaviors in these two examples, which have 
different implications in terms of dark energy. Imposing the vanishing of 70 
at A ew produces a very fast growth of 70(A) at larger energies A ew < A < 
Anni/i while imposing the vanishing at A uni f produces a 70(A) that quicky 
decreases to a large negative value and then slowly decreases to become 



small again near A eiu , as illustrated in the two examples of Figure 12 A 
negative cosmological constant adds to the attractive nature of gravity hence 
it counteracts other possible inflation mechanism that may be present in 
the model, while a positive cosmological constant counteracts the attractive 
nature of gravity and can determine expansion. Thus, we see in these two 
examples that different choices of the free parameter in the model can lead 
to very different qualitative behaviors of the effective cosmological constant. 

The range of different variable cosmological constants can be further ex- 
panded by considering other possible values of the parameter f± in the sur- 
face 7o(A, fi, $2) associated to a chosen value of fi- 

Again, as we have seen in the case of the effective gravitational constant, 
we have here two distinct mechanisms for vacuum-decay: the running 70(A) 
as described above, depending on the parameters j% and and on the 
renormalization group equations, and a further modification due to the non- 
minimal coupling with the Higgs field, as in |25j . 

Namely, if one has an unstable and a stable equilibrium for \H\ 2 , one 
can encounter for the cosmological constant the same type of phenomenon 
described for the effective gravitational constant by the gravity balls, where 
the value of 70,^ is equal to 70 at the unstable equilibrium \H\ 2 = on some 
large region, while outside that region it decays to the stable equilibrium 
\H\ 2 = v 2 for which 



7o(A) 
1 - 167rG eff (A)eo^ 



7o,tf(A) = — _ u ._ r< , v w 



This behaves differently in the case of our system from the model of [25] . In 
fact, here the non-minimal coupling is always the conformal one £0 = 1/12, 
while the gravitational constant is itself running, so we obtain the following. 
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Figure 13. Running of 7o,/r(A) with f 2 
7o(A uni /) = 0. 



1 and with 



Proposition 4.13. The non-minimal conformal coupling of the Higgs field 
to gravity, 

-Co / R\H\ 2 ^-gd A x 

changes the running of the effective cosmological constant to 

7o (A) 



(4.51) 
This gives 



7o,h(A) 



7o,tf(A) 



-167rG cfi (A)^\H\ 2 ' 
7o (A) 

4a(A)(2/ 2 A2a(A)-/ e(A)) 
A(A)b(A)(192/ 2 A2-2/ c(A)) 



for |iJ| 2 ~/x§/(2A ). 



Proof. Arguing as in [25J we obtain (4.51 ), which now depends on the run- 



ning of both 70(A) and G e g(A), as well as on the Higgs field \H\ 2 . Assuming 
the latter to be nearly constant with \H\ 2 ~ /Xq/(2Aq), we obtain 

GcsM A ) _ 7o (A) 



(4.52) 7o>H (A) = 70(A)- 



Geff(A) 



1 



3 



r «m l 2/2A2a(A)-/oc(A)q(A) 
^eff^J [ ^A(A)b(A) 



where, as in the treatment of the variable effective gravitational constant, 
we have 

■ 2 2/ 2 A 2 a(A)-/ e(A)a(A) 



\H\ 



2A 



vr 2 A(A)b(A) 



using the ansatz ( |4.5[ ) for A (A), and G e ^ H (A) = G eff (A) 



□ 



To see how the interaction with the Higgs field can modify the running of 
the effective cosmological constant, consider again an example that exhibits 



the same behavior as the second graph of Figure 12 This is obtained, 
for instance, by choosing / 2 = 1 and f^ satisfying the vanishing condition 
7o(A nri j/) = 0. Then the running of 7o,#(A) given in (4.52) exhibits an 



overall behavior similar to that of 70(A) but with a change of sign and a 



singularity, as shown in Figure 13 
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Figure 14. Running of 7 (A)/(87rG e ff (A)) with / 2 = 1 and 
with 7o(A TOi /) = 0. 



4.13. Early time bounds estimate. In models with variable cosmological 
(and gravitational) constant, some strong constraints exist from "early time 
bounds" on the vacuum-energy density that are needed in order to allow 
nucleosynthesis and structure formation, see [7J, [27] . 

In particular, in our model, this means that we can prefer choices of the 
parameters / 2 and for which the ratio 

(a 7o (A) 

(453) ^GW(A) 

is small at the electroweak end of the energy range we are considering, so 
that the resulting early time bound will allow the standard theory of big- 
bang nucleosynthesis to take place according to the constraints of [7j , [27] . 
Notice that, since we have 7o,iy(A)/G e ff j #(A) = 7o(A)/G e g(A) the esti- 



mate on (4.53 ) is independent of further effects of interaction with the Higgs 
field, and it only depends on the choice of the parameters / 2 and of the 
model. 

This type of estimate can be used to select regions of the space of param- 



eters / 2 and f± that are excluded by producing too large a value of (4.53) 
at A = A ew . For example, consider the two cases of Figure [12] where we 
set / 2 = 1 and we choose so that it gives vanishing 70(A) at A ew or at 
A-unif- The first case gives g T'wT' 1 ^ = 0, while the second gives a very 

■I 07T(_r c ff' \l\-ew ) 

large negative value of (|4.53|) at A ew of -5.93668 x 10 86 with a running of 



(4.53) as in Figure 14 



5. Dark matter 

We have concentrated in this paper on early universe models, with an em- 
phasis on various mechanisms for inflation and dark energy that arise within 
the NCG model, and their mutual interactions. Another main question 
about cosmological implications of the noncommutative geometry models of 
particle physics is whether they can accommodate possible models for cold 
dark matter. This exists the domain of validity of the asymptotic expan- 
sion of the spectral action, as such models apply to a more modern universe 
than what is covered by the perturbative analysis of the spectral action 
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functional. However, one can at least comment qualitatively on possible 
candidates within the NCG model for dark matter particles. In addition to 
the minimal Standard Model, one recovers form the computation of the as- 
ymptotic formula for the spectral action and from the fermionic part of the 
action in [17 additional right handed neutrinos with lepton mixing matrix 
and Majorana mass terms. Thus, since this is at present the only additional 
particle content beyond the minimal standard model that can be accom- 
modated in the NCG setting, it is natural to try to connect this model to 
existing dark matter models based on Majorana mass terms for right handed 
neutrinos. 

Those that seem more closely related to what one can get within the NCG 
model are the ones described by Shaposhnikov-Tkachev |45j . Shaposhnikov 
|44j . and Kusenko [33] • In these dark matter models, one has the usual 
active neutrinos, with very small masses, and an additional number of sterile 
neutrinos with Majorana masses. In the case of the z^MSM model of |44j . 
one has three active and three sterile neutrinos. In these models the sterile 
neutrinos provide candidate dark matter particles. However, for them to 
give rise to plausible dark matter models, one needs at least one (or more) 
of the sterile neutrino Majorana masses to be below the electroweak scale. 
In the detailed discussion given in [33] one sees that, for example, one could 
have two of the the three Majorana masses that remain very large, well above 
the electroweak scale, possibly close to unification scale, while a third one 
lowers below the electroweak scale, so that the very large Majorana masses 
still account for the see-saw mechanism, while the smaller one provides a 
candidate dark matter particle. 

It is possible to obtain a scenario of this kind within the NCG model, pro- 
vided that one modifies the boundary conditions of [1] in such a way that, 
instead of having three see-saw scales within the unification and the elec- 
troweak scale, with the smallest one already at at very high energy around 
10 12 GeV, one sets things so that the lowest Majorana mass descends below 
the electroweak scale. A more detailed analysis of such models will be car- 
ried out in forthcoming work where we analyze different choices of boundary 
conditions for the RGE flow of the model. 

6. Conclusions and perspectives 

We have shown in this paper how various cosmological models arise nat- 
urally from the asymptotic expansion of the spectral action functional in 
the noncommutative geometry model of particle physics of [T7] and the run- 
ning of the coefficients of this asymptotic expansion via the renormalization 
group equations of [T] . 

We have seen in particular the spontaneous emergence of conformal grav- 
ity and Hoyle-Narlikar cosmologies at phase transitions caused by the run- 
ning of the effective gravitational constant. We described effects of this run- 
ning on the gravitational waves and on primordial black holes. We described 
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mechanisms by which this running, combined with the conformal coupling 
of gravity to the Higgs field, can generate regions of negative gravity in the 
early universe. We discussed the running of the effective cosmological con- 
stant and slow-roll inflation models induced by the coupling of the Higgs to 
gravity. We discussed briefly the connection to dark matter models based 
on right handed neutrinos and Majorana mass terms. 

A planned continuation of this investigation will cover the following topics: 

• Varying boundary conditions for the RGE flow. 

• Dark matter models based on Majorana sterile neutrinos. 

• Extensions of the NCG model with dilaton field. 

• Nonperturbative effects in the spectral action. 

• Exclusion curve from spectral index and tensor to scalar ratio. 

7. Appendix: boundary conditions for the RGE equations 

We recall here the boundary conditions for the RGE flow of [1] and we 
discuss the compatibility with the condition assumed in [17!. A more de- 
tailed analysis of the RGE flow of [1| with different boundary conditions and 
its effect on the gravitational and cosmological terms will be the focus of a 
followup investigation. 

7.1. The default boundary conditions. The boundary conditions for 
the RGE flow equations we used in this paper are the default boundary 
conditions assumed in [Tj. These are as follows. 



yu = 0.0000482105 - 3.382 x 10~ 15 i 

yi2 = 0.000104035 + 2.55017 x 10~ 7 i 

yi 3 = 0.0000556766 + 6.72508 x 10~ 6 i 

y 2 i = 0.000104035 - 2.55017 x 10~ 7 i 

y 22 = 0.0005 092 79 + 3.38205 x 10~ 15 i 

y 2 3 = 0.00066992 - 4.91159 x lO" 8 ? 

y 31 = 0.000048644 - 5.87562 x 10~ 6 i 

y 32 = 0.0005 85302 + 4.29122 x 10~ 8 i 

y 33 = 0.0159991 - 4.21364 x lQ- 20 i 




1 



2 




For Y d (A unif ) 



(yij) they have 
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1 

Ynu(Aunif) = | 0.5 

0.1Y 

-6.01345 x 10 14 3.17771 x 10 12 -6.35541 x 10 11 

M(Knif) = ( 3.17771 x 10 12 -1.16045 x 10 14 5.99027 x 10 12 

-6.35541 x 10 11 5.99027 x 10 12 -4.6418 x 10 12 

7.2. Constraints from NCG. There are constraints on the boundary con- 
ditions at unification in the noncommutative geometry model. Those were 
described in |17j and we report them here below. 

• A constraint on the value at unification of the parameter A: 

UA ^ _ ^ KKnif) 

2foa{A unif ) z 

• A relation between the parameter a and the Higgs vacuum: 

Vafo = 2M W 
tt g 

• A constraint on the coefficient c at unification, coming from the see- 
saw mechanism for the right handed neutrinos: 

7 S C(,A nni jJ < 

JO JO 

• The mass relation at unification: 

(7-1) Yl ( m » + m e+ 3m l + 3m d)U=A™/ = 8M^\ A=Aunif , 

generations 

where m u , m e , m u , and are the masses of the leptons and quarks, 
that is, the eigenvectors of the matrices tfji, 6±i, 5^ and 5^, respec- 
tively, and M\y is the W-boson mass. 

Clearly, not all of these constraints are compatible with the default boundary 
conditions of [T]. So either one relaxes some of these conditions, as we have 
been doing in the present paper, or one performs a wider search for more 
appropriate and fine tuned boundary conditions for the RGE flow, analyzing 
how different choices of boundary values affect the behavior analyzed in this 
paper. This is presently under investigation. 
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